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Abstract 

This paper deals with the topology of the relative trajectories in flight formations. The purpose is 
to study the different types of relative trajectories, their degrees of freedom, and to give an adapted 
parameterization. The paper also deals with the research of local circular motions. Even if they 
exist only when the reference orbit is circular, we extrapolate initial conditions to the eccentric 
reference orbit case. This alternative approach is complementary with traditional approaches in 
terms of cartesian coordiantes or differences of orbital elements. 


1 Introduction 

Flight formations have become a key technology for present and future space missions. From a dy- 
namical point of view, they have been studied since the first spatial rendezvous operations [2]. There 
is a certain number of publications in the literature presenting different analytical theories for tempo- 
ral evolution of the relative motion [4], [5]. Different authors have explored the possibility of finding 
periodic relative orbits under the effect of perturbations [1], [7]. We consider that, before taking into 
account perturbations, there is a field of investigation which consists in describing all the possible 
relative trajectories for the non-perturbed case. This paper tackles this task. We lay stress on circular 
or near circular local trajectories because circular local trajectories have special properties which make 
them specially useful for space missions (interferometry missions, LISA mission, ....). 

First, we study the circular reference orbit case. Since the relative trajectory is, under some assump- 
tions, an ellipse, we propose to study it using "local orbital elements" defined in [3]. We are explaining 
that the six local orbital elements are not independent, and that the relative ellipse has only four 
degrees of freedom. Second, when the reference orbit is slowly eccentric, the relative motion is no more 
an ellipse. Our local orbital elements are always a well-suited parametrization because the resulting 
relative motion is a kind of perturbed ellipse. For highly eccentric references orbits, we propose a 
parameterization based, only, on four parameters, two for the inplane motion, and two for the out-of- 
plane motion. 

The paper has the following structure: second section is dedicated to the statement of the problem. 
Third section deals with the topology of relative motion while section four is devoted to the search of 
circular relative motions. 

2 The equations of the relative motion 

2.1 The statement of the problem 

In the whole text, we use the well-known keplerian elements: the semi-major axis a, the eccentricity 
e, the inclination i, the right ascension of ascending node D, the argument of perigee a;, and the mean 
anomaly M. The motion, studied in an inertial reference frame denoted IJK, is described through 
temporal series of keplerian elements as well as of positions ~t^\ijk and velocities ~v\ijk- We use the 
following notations: 
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EO — {a,e,i,VL,uj,M) 

I IJK 


X \IJK = 


V \IJK 


From now on, we will consider a reference orbit which will be described by its orbital elements or by 
its position and velocity. This reference orbit can be the orbit of one of the satellites of the formation 
or it can correspond to a fictitious point. For simplicity, we will name it the reference satellite indicated 
by the subscript r . Inertial position and velocity can also be projected into the orbital local frame 
(. RTN ) defined by the given reference orbit. Relations between projections into IJK frame and RTN 
frame are quite simple: 


X \IJK 



x | RTN 


where the matrix 1Z(EO) is: 


K(EO) 


cos cos u — sin ft sin u cos i 
sin Q cos u + cos Q sin u cos i 
sin u sin i 


— cos sin u — sin Q cos u cos i 

— sin Q sin u + cos Q cos u cos i 

cos u sin i 


sin Q sin i 
— cos Q sin i 
cos i 


with u = cj + / , where / is the true anomaly. 

The relative motion between a satellite and the reference satellite can be described by the difference 
of position and velocity projected in the orbital frame of the reference orbit: 


A X | RTN = X | RTN ~ X t \rtn = 


A r I RTN 
AT*| RTN 


or by the differences of orbitals elements between the two orbits: 


( 1 ) 


A EO = EO — EO r (2) 

We will also use the following notation: 

AX\ rtn = (AR,AT,AN) t 

AV | rtn = (AVr, AVt, AVn) T 


2.2 The temporal evolution of the relative motion 

In order to obtain easy analytical expressions of the relative motion, we study the linear, non-perturbed 
problem. As it is proved in [3], the equations of the motion of this problem are: 


Y CL €!■ 3 TL (3 

A R(t) = — Aa 0 — a r cos / r Ae 0 H — fsin/ r AMo — sin f r (t — to)Aac> (3) 

a r Tj r 2 T) r 

( J ^ \ CL^Tl 3 CL TL T1 

1 H — r— ) sin/ r Aeo + r r cosi r Afio + r r Acjo H — ^AMq — - r r r (t — to)Aao 
T)* a r J r r 2 r r 

A N(t) — r r sin u r Ai$ — r r sin i r cos u r AQo 

These equations give temporal evolution of relative position projected in the local orbital frame 
as function of six parameters: AEO(to). They are equivalent to the Lawden’s equations [6]. Law- 
den’s equations are parameterized with initial relative position and velocity while our equations are 
parameterized with initial differences of orbital elements. Our parametrization gives more compact ex- 
pressions and enables a better physical understanding of the role of each initial condition. Thereafter, 
we describe resulting relative motion in three cases: (i) when the reference orbit is circular, (ii) when 
the reference orbit is slightly eccentric, (iii) when the reference orbit is strongly eccentric. 
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3 The topology of the relative motion 


3.1 The circular reference orbit case 

When the reference orbit is circular, we parameterize the two orbits with a set of non-singular elements: 
ENS — (a, C, i, fi, S, A) T , where the elements (e,cj,M) have been replaced by the elements (C, S, A) 
defined by: 

C — ecoscj (4) 

S — esin u (5) 

A = lj + M ( 6 ) 

We also define the differences of non-singulars elements as: 

A ENS = ENS - ENS r (7) 

The particularization of equations (3), after the introduction of non-singular elements, gives: 


AR(t) = 

Aao — (cos AACo + sin AAS'o) 


( 8 ) 

A T(t) = 

/ 3 

a r ( 2 sin AACo — 2 cos AAS'o + AAo — cos i r AQq — 

- t 0 )Aa 0 ^ 


AN(t) = 

a r (sin AAio — sin i r cos XAQq) 




This solution is equivalent to the well-known Clohessy-Wilthsire equations for the non-perturbed 
motion. We note that these equations correspond with the standard parametrical equation of an ellipse 
centered on the origin of the local frame except for two kind of terms: 

• terms providing from a difference of semi-major axis : a different semi-major axis gives different 
orbital frequencies which produce a secular growth of the difference in the T component. As this 
effect destroys the formation in a short period of time, we do the hypothesis: Aao = 0. This 
hypothesis may not be true for docking or rendezvous operations, but it is always fulfilled in 
flight formations. 

• constant term on the T axis : these terms, (AAo — cosi r Afio ) 5 can be removed just by a shift on 
the origin of the axis. We do not take care of these terms. 

Disgarding foregoing terms, equations ( 8 ) rewrite: 

A R — — a r ACocos A — a r ASosin A (9) 

AT = — 2a r ASocos A + 2a r ACosin A 
AN = — a r AJ4ocos A + a r Aiosin A 

which correspond to the elliptical motion. As it is usually done on the two-body problem, the 
elliptical motion will be parameterized through a set of orbital elements called local orbital elements 
( eoi ): semi-major axis (a/), eccentricity (e/), inclination (i/), longitude of ascending node (O/), longitude 
of perigee (cj[) and the phase (cpi), defined as the difference of phase between the reference anomaly, 
and the anomaly of the local orbit. 

There are two main differences with respect to classical keplerian elliptical motion (i) the origin of 
the axis does not correspond to a focus of the ellipse, but with the center. This leads to an ambiguity 
on the definition of the longitude of the perigee which is solved establishing that u E [ 0 , tt] (ii) the 
angular velocity is not dependent on the distance to the origin, it is constant and equal to the angular 
velocity of the reference frame. That is, the period of any local orbit corresponds to the period of the 
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Local orbital elements 

a, : semi -maj or axi s 

e, = 1 - (b/a,) 2 : eccentri ci ty 

i, : inclination 

Q, : longitude of ascending 
node 

C 0 | : longi tude of perigee 
M, : anomaly 


Figure 1: definition of local orbital elements 


reference frame. 

It is possible to obtain relations between the local orbital elements and the initial conditions as it is 
done in [3]. It is also possible to prove that only four local orbital elements are independents. The size 
(i ai ) and the constant phase (■ ipi ) can always be chosen to our convenience. The other four parameters 
may be separated in two groups: the elements which give the form of the local orbit (e/,a;/), and the 
elements which give the orientation of the plane of the local orbit (i/,0/). One group determines the 
other. We have decided to express the local eccentricity and the local perigee as function of other 
variables (i/,0/). Tedious but simple algebraic manipulations yield: 

ef [9 + 6 tan 2 i/(4cos 2 0/ — sin 2 O/) + tan 4 i/(4cos 2 0/ + sin 2 O/) 2 ] 

(2 — e ?) 2 (5 + tan 2 i/(4cos 2 0/ + sin 2 O/)) 2 

and: 


o 3(cos 2 0/ 
tan 2 uji 


sin 2 O/) + sin 2 i/(cos 2 0/ + 4 sin 2 O/) 
3 cos i/ cos 0/ sin 0/ 


tana;/ — 1 = 0 


(ii) 


3.2 The low eccentric reference orbit case 

When the reference orbit is sligthly eccentric, we simplify equations (3) by means of an expansion in 
powers of the eccentricity up to first order: 


A R 

CLy 

AT 

a r 


- cos A (AC 0 + SVAA 0 ) - sin A (A S 0 - C r AA 0 ) + 0(e 2 ) (12) 

sin A [2ACo + S r (AAo — AOq cos i r )\ + cos A [— 2ASo + C r ( AAo — A Oq cos i r )\ 

+ (AAo — AOo cosi r ) — Aeo^- (sin2Acos2a; r — cos2Asin2a; r ) + 0(e 2 ) 
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Figure 2: Eccentricity as function of the inclination and the longitude of the ascending node 
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Figure 3: Perigee as function of the inclination and the longitude of the ascending node 
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= sin AAffi — cos A sin i r AQq + ^ (cos c o sin i r AQq — sin u r Aio) 

Qjy> Z 

g g 

+ -^ cos2A (cos ay sini r A^o + sinoyAffi) + sin2A (sin ay sini r Afio — cosoyAffi) + 0(e 2 ) 


In this set of equations we identify following terms: 


• constant terms : there are constant terms not only in T axis (as in the circular case) but also in 
the N axis. Once again, they can be cancelled by changing the origin of the axis. 

• elliptical terms : the terms in cos A and sin A. The coefficients differs from the coefficients of the 
circular reference orbit case. The local orbital elements can be computed using the same relations 
as in the circular reference orbit case but with new parameters P : 


A — — a r (ACq "T SVAAo) B — —a r ( ASq — C r AAo) 

C — a r (—2 ASq + C r (AAo — Afio cos i r )) D — a r ( 2ACo + S r (AAo — cos i r )) 

E — — a r sinz r AJ2o F — a r AiQ 

• double orbital frequency terms : relative motion is no more an ellipse because of these terms which 
do not correspond to an elliptical motion. 

From a topological point of view there is a major modification with respect to the circular case: 
A A plays a new role. While in the circular case it determines only a constant, in the low eccentricity 
case it plays a role on the ellipse determination and on the double orbital frequency terms. The 
modifications with respect to circular case are of the order of the eccentricity, since the eccentricity is 
low, the difference are not very important. 

3.3 The high eccentric reference orbit case 

Our departure point is equations (3). For the same precedent reasons, we impose A a — 0, obtaining 
equations: 


U/v. Cj’ 


A R(t) = — a r cos f r Aeo H sin f r AMq 


T)r 


(13) 


1 ry 


AT (t) — a r 1 H — 77 — sin f r Aeo + r r cos i r AJ2q + r r Acjq + 

' rtf a r 1 


a^rj r 


AM n 


A N(t) = r r sin u r Azq — r r sin i r cos u r AQq 

The in-plane and out-of-plane motions can be decoupled doing the following changing of variables: 


A uJ — Au + AQ cos i r (14) 

AQ' = A ft sini r 

After, we break down the variable u = uj + / into his components: 


A i" = sino; r Ai — coscj r Ai2 / (15) 

A ft" — coscj r Ai + smu; r Afl' 


Finally, we obtain: 
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X 


cos f r + eK\ sin f 7 


(16) 


AR 

a r Ae 
AT 
a r Ae 
AN 
a r Ae 

with the coefficients: 


2 + e cos f r . , , . v . 

— — sin / r + Ki (1 + e cos f r ) + 

1 + e COS J r 

— — — (K 3 cos fr + Ki sin f r ) 

1 + e COS Jr 


K 2 

1 + e cos f r 


Kf 

K 2 

K a 


AM 

Vr 

r ] 2 ( Auo + cos i r AQ) 

rft (sinoyAi — coso; r sinz r Afl) 

rft (coscj r Ai + sincj r sinz r Afi) 


(17) 


These equations express all the possible relative trajectories as function of four parameters. K\^K 2 
parameterize the in-plane motion and the out-of-plane motion. The out-of-plane motion is 

given by the sum of a sinus and a cosinus functions divided by (1 + ecos f r ). The in-plane motion is 
more complicated. K\ produces a circular motion: 


x — — cos f r + eKi sin f r 

y — K\ — eK\ cos f r + sin f r 


(18) 


and K 2 produces a motion only on y axis: 


y ~ K 2 


—eK 2 cos f r 
1 + e cos f r 


The values of these constants determine the form of the in-plane motion. In figure (4) we plot the 
form of the in-plane motion for different values of K\^ K 2) when the reference orbit eccentricity is 0.6. 
We verify that for high values of K\ the motion is circular while K 2 gives linear motions on T axis. 


4 Looking for local circular motions 

The planes where the local motion is circular, have a very interesting property: due to the fact that 
the local orbital frequency is independent of the distance to the origin, all the satellites rotate at 
the same angular velocity around the center of the circle and their relatives distances do not change. 
As a consequence, all initial configurations remain invariant all along the trajectory. These planes 
represent very interesting possibilities for flight formation where it is necessary to keep inter-relative 
distances constant, as in interferometry missions or the future LISA mission. Actually, residual varia- 
tions of the distance will be produced by perturbations and non-linear effects that have been neglected. 


4.1 When the reference orbit is circular 

There are two points in figure (2) where local eccentricity is zero. Using equations (10) and (11), it is 
possible to determine their local plane: 


111 

112 
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60° 

-60' 
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fin = 90° 
fl n = 90' 




Figure 4: In-plane motion as function of K\ and K 2 


with the corresponding initial conditions: 


A ENSi = (0, AC, -VSAS, VSAC, AS, AA) T 

> (19) 

AENS 2 = (0, AC, VSA S, — Vs AC, AS, AA) T 

We recall that, due to equations (9), AC and AS determine the size (a/) and the constant phase 
of the ellipse (cpi). 


4.2 When the reference orbit is slightly eccentric 

When the reference orbit is not circular, the set of equations (12) reveals that it is not possible to find 
perfect circular local motions. But, our analysis of these equations enables to establish the necessary 
conditions to provide local motions as close as possible to a circle. When changing the parameterization 
(12) by using differential orbital elements we obtain: 

= — cos/ r Ae + e r sin / r AM (20) 

= 2 sin f r Ae — e r cos f r sin f r Ae + (AM + A a/) + e r cos f r (AM — A a/) 

= (A i r cos f r + AfT sin f r ) — e r cos f r (A i r cos f r + AQ' sin f r ) 
e introduced the following variables: 


A R 


AT 

CLr 

AN 

a r 


where we hav 


Ao/ — Acj + Aflcosi r (21) 

A i — sinuj r Ai — cosuj r smi r AQ 
An' = cosuj r Ai + sino; r sin z r All 

Considering only elliptical terms, the following initial conditions enable to obtain the same circular 
motion as in the circular reference orbit case (same local inclination and ascending node): 

• SET 1 


A a — 0 

a ^ x /77 COS LOj- . 

An = -V3 -Ae 

sin i r 


Ae = Ae 


A lj = v/3 


COS LJ r 

tan i r 


Ae 


A i = V3sm(jj r Ae 
AM = 0 


• SET 2 


A a = 0 

AQ = VSeTr—^-AM 

sin i r 


Ae = 0 
A uj — — ( 1 v/3' 


sin (jj r 
tan i r 


AM 


A i = V3e r cos (ju r AM 
AM = AM 


The first set of initial conditions is parameterized by a difference of eccentricity, while the second 
set is parameterized by a difference of anomaly. Thereafter, we refer to the set 1 as "difference of 
eccentricity" and to set 1 as "difference of anomaly". Resulting motion is not circular because of 
double orbital frequency terms. In figure (5) we compare the variation of the distance from the origin 
for each set of initial conditions. Both sets of conditions produce double orbital frequency terms in the 
N axis, but, Ae produces also these terms in the T axis. That’s why the variations are more important 
in the left part of the figure (5). 
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Distance to the origin 



Distance to the origin 



Figure 5: Distance to the origin for different initial conditions (left: difference of eccentricity , right: 
difference of anomaly) 


4.3 When the reference orbit is highly eccentric 

It is possible to arrange equations (13) to obtain following form: 


A R 

dj* 

AT 

eLy 

AN 

a r 


— cos / Ae H — sin f AM 

T] 

( AM 

2 sin f Ae + e cos / I 

Tj 2 


(22) 


tj 2 Auj' 


1 + e cos / 


r] 1 + e cos / 
(A i cos / + AD" sin /) 


+ ( Tj 2 Aj + 


AM 


cos / sin / 

e: Tb r Ae 

1 + e cos f 


This arrangement allows the comparison with the low eccentric form of the equations 20. Comparing 
the two sets, we identify the transformation of elliptical terms when the eccentricity of the reference 
orbit grows. We identify the following terms: 


— cos / + e sin / 

2 sin f Ae + e cos / (AM — Ac</) 
(A i! cos / + A ft" sin /) 


— cos / H — sin / 
rj 

n . „ . (AM 

2 sin / Ae + ( 

V V 

T] 2 (A i! cos / + AD" sin /) 


— r] 2 A(jj' ) ecos / 




We impose the same previous conditions to new terms. Resulting motions are not circular because 
of double orbital frequency terms and non linear terms in eccentricity, that become very important. 
Anyway, these initial conditions produce a class of relative motions that is as close as possible to a 
circle. 


• SET 1 


A a = 0 

. ^ y/3 COS L O . 

AD = — Ae 

rj z sin i 


Ae = Ae 

. v^3 cos uj A 

Acj = -V Ae 

rj z tan i 


A • \/3 . . 

A i — sinwAe 

AM = 0 


• SET 2 

A a = 0 

AQ = VS 


Ae = 0 


e r sin uo r 
r] 3 sini r 


AM 


. 1 / . nz sin a; 

Acj = o 1 + v3e r 7 


tan L 


A i — a/3— ^ cosu; r AM 

Tj 6 


AM AM AM 
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In the following figures (6) and (7) we have drawn the evolution of the circular motion when the 
eccentricity of the reference orbit grows. The other parameters of the reference orbit play no role on 
the form of the trajectory (for the keplerian case). Figures (6) and (7) show how, for very big reference 
orbit eccentricities, motion is far from being circular. 

Local Relative Motion 


e = 0 
e = 0.1 
e = 0.3 



Figure 6: Evolution of the local circular motion produced by a difference of eccentricity when the 
eccentricity of the reference orbit grows 


Local Relative Motion 


e = 0 
e = 0.1 
e = 0.3 



Figure 7: Evolution of the local circular motion produced by a difference of anomaly when the eccen- 
tricity of the reference orbit grows 
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5 Conclusions 


This paper explores a field in relative motions which is nowadays not well-known . We characterize all 
the possible relative motions for all type of eccentricity of the reference orbit. We minimize the number 
of relevant parameters thanks to a detailed analytical study. In the case of the circular reference orbit 
case, we take the advantage of the representation in terms of local orbital elements. In the case of 
very high eccentric reference orbit it is not possible to use the local orbital elements anymore, but we 
separate the motion as a sum of four elementary effects. 

Local circular motion has very interesting properties for space purposes. That’s why we have identified 
the configurations that produce this motion in the circular reference orbit case. When the reference 
orbit is eccentric, there are no more local circular motions, but we have found the motions that are 
closest to a circle. 

Future work will be centered in the introduction of perturbations and their effects on the relative 
motion. Another interesting work is to find the correct parametrization of relative motion when 
reference orbit is highly eccentric. 
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